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RESUMEN 

Se presenta el analisis sfsmico, mediante el Metodo Espectral de los 
bloques estructurales 6 y 7 de la Universidad de Fuerzas Armadas ESPE, que 
tiene aisladores de triple pendulo de friccion colocados sobre los plintos de 
cimentacion y sobre las columnas del primer piso. 

Se realiza un analisis espacial, considerando tres grados de libertad por 
planta: dos componentes de desplazamiento horizontal y una rotation de piso con 
respecto a un eje perpendicular a la losa. Pero para mostrar el uso del sistema de 
computation CEINCI-LAB tambien se realiza un analisis piano. 

El analisis se lo realiza tanto en sentido X, como en sentido Y, para dos 
espectros: el de diseno DBE y el maximo considerado MCE. La rigidez efectiva del 
aislador FPT8833/1 2-1 2/8-6 se la obtuvo mediante el modelo de Me Vitty y 
Constantinou (2015) trabajando con las dos fases para el sismo DBE y con las tres 
fases para el sismo MCE. 

Se comparan los desplazamientos laterales hallados en el analisis 
espacial, con los encontrados en el analisis piano, tanto para un portico en sentido 
X, como para un portico en sentido Y, con el sismo DBE. 

Finalmente, se comparan los resultados maximos probables de la 
estructura aislada con los que se hallan al no considerar aisladores sismicos y se 
observa como se mejora notablemente el comportamiento al tener la estructura 
aisladores FPT. 

Palabras Claves.- Aisladores de triple pendulo de friccion. Metodo Espectral. 
Sismo de Diseno. Sismo Maximo Considerado. 
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ABSTRACT 

The seismic analysis using the modal spectral method of the structural 
blocks 6 and 7 of the Armed Forces - ESPE University is presented, This structural 
blocks have triple friction pendulum isolators placed on foundation plinths and over 
the columns on the first floor 

A spatial analysis is performed, considering three degrees of freedom 
per each floor: two components of horizontal displacement and a rotation of the 
floor with respect to an axis that is perpendicular to the slab. But to show the use 
of the computer system CEINCI LAB also a flat analysis is performed. 

The analysis is performed both in the X and Y direction, for two 
spectra: the design DBE and the maximum considered MCE. The effective 
stiffness of the isolator FPT8833/1 2-1 2/8-6, was obtained by the model of Me Vitty 
and Constantinou (2015) working with the two phases for the DBE earthquake and 
the three phases for MCE earthquake. 

The lateral displacements found in spatial analysis are compared with 
those found in the flat analysis. For a frame in X direction, as a frame in Y direction, 
considering the DBE earthquake. 

Finally the maximum probable results of the isolated structure are 
compared with those who are not considering seismic isolators and it shows how 
the behavior is markedly improved as the structure have the FPT isolators. 

Keyword. - Triple friction pendulum isolators. Spectral method. Design 
Earthquake. Maximum Considered Earthquake. 

1. INTRODUCTION 

En el Proyecto Arquitectonico del nuevo campo de investigation cientffica 
y de post grados de la Universidad de Fuerzas Armadas ESPE, se tienen dos 
bloques estructurales identificados como 6 y 7. El primero de ellos esta destinado 
al Centro de Investigaciones de Electronica y Seguridad y Defensa; en cambio en 
el bloque 7 se tienen las aulas de post grado y laboratories de informatica. 

Una vez que se decidio utilizar aisladores sfsmicos, estos dos bloques 
estructurales desde el punto de vista arquitectonico, fueron unidos para evitar la 
junta de construction de aproximadamente 1.20 m, que conlleva tener dos 
bloques estructurales con aisladores sfsmicos. Al tener un solo bloque se elimina 
la junta y con ello se aporta a mejorar las fachadas del proyecto arquitectonico, 
debido a que se tiene una estructura continua. 

En la figura 1 se presenta, las plantas que conforman este bloque 
estructural que se lo va a identificar como 67, ya que es un solo bloque. En la 
parte superior izquierda se tiene el subsuelo, a la derecha de la losa de aislacion; 
posteriormente en la parte inferior se observa a la izquierda las losas de los 
niveles 4.94 y 9.88; al centra del nivel 14.82 y a la derecha del nivel 19.76; todos 
estos niveles tienen un voladizo a la derecha. 

En la figura 2 se presenta un code longitudinal del bloque 67, en que se 
aprecian los aisladores sobre la cimentacion en los tres primeros ejes y sobre las 
columnas del subsuelo en los siguientes ejes. 
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Figura 1 Plantas del Bloque Estructural 67 de la UFA-ESPE. 
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Figura 2 Corte longitudinal del Bloque Estructural 67. 



2. MODELO DE TRES FASES DEL AISLADOR FTP 

Para los aisladores de triple pendulo de friccion, Fenz y Constantinou 
(2007, 2008), Fadi y Constantinou (2010) desarrollaron un modelo de 5 fases para 
definir el comportamiento no lineal de un aislador FPT; considerando para ello 12 
parametros geometricos y 4 coeficientes de race para el aislador. 

Normal mente se construyen aisladores de tal manera que de los 12 
parametros se reducen a 6 y de los 4 coeficientes se bajan a 2. Para estas 
condiciones. Me Vitty y Constantinou (2015) desarrollaron un modelo de 3 fases, 
pero antes de ello se indica a continuation los 6 parametros geometricos. 

• = i?4 ^ f?2 = 

• y d 2 = d 2 

• = h 4 y h 2 = h 3 

El significado de las variables se indica en la figura 3, en que se presenta 
la geometrla del aislador FPT8833/1 2-1 2/8-6, que se empleo en el Proyecto. 

Los dos coeficientes de roce son: n ± = /i 4 y /i 2 =A 3 - Para estas 
condiciones se tienen 6 parametros geometricos y 2 condiciones de roce, con lo 
que el modelo de 5 fases se reduce a un modelo de 3 fases, las mismas que se 
indican en la figura 4; en la parte superior se presenta la fase II y en la parte 
inferior la fase III. 

Una breve descripcion de las fases del modelo, es la siguiente; en la 
primera fase el deslizamiento se produce en el aislador interior de radio R 2 que es 
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igual afi 3 y que tiene un coeficiente de race ju 2 . El desplazamiento de esta fase es 
muy bajo. 

i 36* , 




Figura 3 Aislador FPT8833/1 2-1 2/8-6, utilizado en el Proyecto 
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Figura 4 Fases II y III del modelo de Me Vitty y Constantinou 2015. 

A la izquierda de la figura 4 se indica el comportamiento del aislador para 
la segunda fase, se aprecia que el aislador interior se desliza sobre las placas de 
radio R t que es igual al R 4 y que tiene coeficiente de race n v En la fase III, el 
aislador interior se topa con los bordes, experimentando grandes desplazamientos. 
Para el sismo de diseno DBE la rigidez efectiva k ef2 se halla con la fase II, y para 
el sismo MCE se encuentra k ef3 con el modelo de fase III pero trabajando con el 
mismo valor de la fuerza horizontal de la fase II con el desplazamiento vertical de 
la fase III. Aguiar y Morales (201 6). 




Desplazamiento 




Desplazamiento 
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3. CARGAS Y RIGIDECES DE LOS AISLADORES 



/c Cl 

La carga permanente D para los pisos interiores es de D = 800 y para 
la cubierta D = 600 Para el analisis slsmico se trabajo con D + 0.25 L. En la 

figura 5 se presenta el mosaico de cargas con las que se hallo la carga vertical 
que Mega a los aisladores, en base de la cual se determinan las rigideces y 
factores de amortiguamiento del modelo de 3 fases. 
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Figura 5 Mosaico de cargas para encontrar la carga que Mega a cada aislador 
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En las tablas 1 a la 9 se indica la carga total que Mega a cada aislador, las 
fuerzas y desplazamientos de la fase II y III del modelo de Me Vitty y Constantinou 
(2015) y la rigidez efectiva para el sismo DBE y para el sismo MCE. Un dato que 
no se ha indicado todavfa es el coeficiente de roce, se calculo con ju 2 = ju 3 = 0.01; 
y con Hi = A 4 = 0.07 



Tabla 1 Calculo de la rigidez efectiva de los aisladores del Portico 32 



Aislador 


W 

(T.) 


Fmax2 

(T.) 


92 

(m.) 


93 

(m.) 


kef 2 
(T/m) 


kef 3 

(T/m) 


1 


86.09 


16.30 


0.5366 


0.566 


30.37 


28.78 


2 


174.79 


33.09 


61.67 


58.43 


3 


174.79 


33.09 


61.67 


58.43 


4 


140.68 


26.63 


49.63 


47.03 



Tabla 2 Calculo de la rigidez efectiva de los aisladores del Portico 33 



Aislador 


W 

(T.) 


F max2 

(T.) 


92 

(m.) 


93 

(m.) 


k e f 2 
( T/m) 


k e f 3 

(T/m) 


1 


180.67 


34.20 


0.5366 


0.566 


63.74 


60.40 


2 


366.81 


69.44 


129.41 


122.63 


3 


366.81 


69.44 


129.41 


122.63 


4 


295.23 


55.89 


104.16 


98.70 



Tabla 3 Calculo de la rigidez efectiva de los aisladores del Portico 36 



Aislador 


W 

(T.) 


F max2 

(T-) 


92 

(m.) 


93 

(m.) 


k e f 2 
(T/m) 


k e f 3 

(T/m) 


1 


189.15 


35.81 


0.5366 


0.566 


66.73 


63.23 


2 


384.03 


72.70 


135.49 


128.38 


3 


384.03 


72.70 


135.49 


128.38 


4 


316.34 


59.89 


111.61 


105.75 



Tabla 4 Calculo de la rigidez efectiva de los aisladores del Portico 39 



Aislador 


W 

(T.) 


F max2 

(T-) 


92 

(m.) 


93 

(m.) 


k e f 2 
(T/m) 


k e fs 

(T/m) 


1 


187.79 


35.55 


0.5366 


0.566 


66.25 


62.78 


2 


381 .27 


72.18 


134.52 


127.46 


3 


381 .27 


72.18 


134.52 


127.46 


4 


326.11 


61.74 


115.05 


109.02 



El factor de amortiguamiento promedio que se hallo para el sismo DBE es 
f eq = 0.2028 y para el sismo MCE es 0.1962. 
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Tabla 5 Calculo de la rigidez efectiva de los aisladores del Portico 41 



Aislador 


W 

(T.) 


Fmax2 

(T.) 


92 

(m.) 


93 

(m.) 


kef 2 
(T/m) 


k e f 3 

(T/m) 


1 


163.18 


30.89 


0.5366 


0.566 


57.57 


54.55 


2 


331.30 


62.72 


116.89 


110.76 


3 


331.30 


62.72 


116.89 


110.76 


4 


294.72 


55.80 


103.98 


98.53 



Tabla 6 Calculo de la rigidez efectiva de los aisladores del Portico 43 



Aislador 


W 

(T.) 


F mcix2 

(T-) 


92 

(m.) 


93 

(m.) 


k e f 2 
(T/m) 


k e f 3 
(T/m) 


1 


69.96 


13.24 


0.5366 


0.566 


24.86 


23.39 


2 


142.05 


26.89 


50.12 


47.49 


3 


142.05 


26.89 


50.12 


47.49 


4 


109.65 


20.76 


38.69 


36.66 


Tabla 7 Calculo de la rigidez efectiva de los aislad 


ores del Portico U 


Aislador 


W 

(T.) 


F max2 

(T.) 


92 

(m.) 


93 

(m.) 


kef 2 

( T/m) 


kef 3 

(T/m) 


1 


86.09 


16.30 


0.5366 


0.566 


30.37 


28.78 


2 


180.67 


34.20 


63.74 


60.40 


3 


189.15 


35.81 


66.73 


63.23 


4 


187.79 


35.55 


66.25 


62.78 


5 


163.18 


30.89 


57.57 


54.55 


6 


69.96 


13.24 


24.68 


23.39 



Tabla 8 Calculo de la rigidez efectiva de los aisladores del Portico V y V’ 



Aislador 


W7 

(T.) 


Fmax 2 
(T.) 


92 

(m.) 


93 

(m.) 


k e f2 

(T/m) 


k e fs 

(T/m) 


1 


174.79 


33.09 


0.5366 


0.566 


61.67 


58.43 


2 


366.81 


69.44 


129.41 


122.63 


3 


384.03 


72.70 


135.49 


128.38 


4 


381.27 


72.18 


134.52 


127.46 


5 


331.30 


62.72 


116.89 


110.76 


6 


142.05 


26.89 


50.12 


47.49 



Tabla 9 Calculo de la rigidez efectiva de los aisladores del Portico W 



Aislador 


W 

(T.) 


Fmax 2 

(T.) 


9 2 

(m.) 


93 

(m.) 


k e f2 

(T/m) 


k e f3 

(T/m) 


1 


140.68 


26.63 


0.5366 


0.566 


49.63 


47.03 


2 


295.23 


55.89 


104.16 


98.70 


3 


316.34 


59.89 


111.61 


105.75 


4 


326.11 


61.74 


115.05 


109.02 


5 


294.72 


55.80 


103.98 


98.53 


6 


109.65 


20.76 


38.69 


36.66 
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4. ESPECTROS DE ANALISIS 

El espectro de diseno DBE, es el que reporta la Norma Ecuatoriana de la 
Construction, NEC-15 para un perfil de suelo tipo C; y el espectro maximo 
considerado MCE se halla multiplicando el espectro de diseno por 1.4 
(Constantinou et al. 2016). Estos espectros se presentan con Ifnea azul en la 
figura 6. 




Ahora bien, se considero un factor de reduction de las fuerzas sfsmicas 
R = 1; un factor de irregularidad en planta 0 P = 1; un factor de irregularidad en 
elevation 0 e = 0.9. Por lo tanto, el espectro elastico para 5% de amortiguamiento, 
en la zona de perfodos menor a 2 segundos se divide para R 0 P 0 e = 0.9; con esto 
en lugar de disminuir el espectro elastico se incrementa como se observa en la 
Ifnea en rojo. 

Para la zona de perfodos largos, mayores a 2 segundos el espectro 
elastico se divide para el valor de B, el mismo que es funcion del factor de 
amortiguamiento de los aisladores. Constantinou etal. (2011). 




Cada aislador tiene un factor de amortiguamiento equivalente. El valor que 
interviene en la ecuacion (1) es el promedio de los factores de amortiguamiento. 
Se encuentra el valor f eq , tanto para el sismo DBE como para el sismo MCE, en 
base a la II y III fase del modelo de Me Vitty y Constantinou (2015). Estos valores 
ya fueron indicados. 

El perfodo de 2 segundos se debe a que la estructura con aisladores 
tienen un perfodo mayor a 3 segundos y se sabe que fundamentalmente el 
sistema de aislacion trabaja en los primeros modos. De tal manera que las 
aceleraciones espectrales para perfodos mayores a 2 segundos son para el 
sistema de aislacion. 
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5. SECCIONES DE LOS ELEMENTOS ESTRUCTURALES 

Se denomina superestructura a los elementos estructurales que se hallan 
sobre los aisladores, y subestructura a los que se encuentran bajo los aisladores. 
Con esta acotacion los elementos de la superestructura son de acero, las 
columnas son tubulares cuadradas de 550/550/20 mm, las vigas son tipo “I” de 
550/300/25/12.5 mm; estas secciones con sus respectivas areas y momentos de 
inercia se indican en la tabla 10. 

Las columnas de la subestructura son de hormigon armado de 80/80 cm. 
En la parte superior de las columnas se tiene una canastilla de 1 .20/1 .20 m, donde 
se asienta el aislador; el hormigon utilizado tiene una resistencia a la compresion 
de 21 0 kg/cm 2 . 

Tabla 10 Secciones de los elementos estructurales del Bloque 67 
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6. MATRIZ DE RIGIDEZ LATERAL DE LOS PORTICOS 

El analisis sismico se realizo utilizando el sistema de computacion 
CEINCI-LAB (Aguiar 2012, 2014). Para ello se numeraron los nudos y elementos 
de cada uno de los porticos, en la forma indicada en las figuras 7 y 8. 

Para la numeracion de los nudos, primero se numeran los restringidos y 
luego todos los demas en cualquier orden. Para los elementos se debe numerar 
primero los de acero; luego los de hormigon y al final los elementos aislador; de 
eta forma se facilita el calculo de la matriz de rigidez total como la contribucion de 
los elementos de acero mas la contribucion de los elementos de hormigon y mas 
la contribucion de los elementos aislador. Si no se desea numerar, los elementos, 
de la forma indicada puede realizarse de otra manera, lo importante es que se 
tengan agrupados los elementos de un mismo material. 

Para ilustrar el procedimiento de calculo de la matriz de rigidez lateral de 
los porticos, en la figura 9 se presenta la numeracion de los grados de libertad del 
portico V. Notese que primero se ha numerado las coordenadas principales 
(laterales) del 1 al 6; luego las coordenadas secundarias: desplazamiento vertical y 
giro. 



Para el portico V se aprecia en la parte superior de la figura 8 que primero 
se ha numerado los elementos de hormigon armado; luego los de acero y 
finalmente los aisladores. 

• Porticos en sentido X 



©000 @ 0 ©© 
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Portico 41 



Portico 43 



Figura 7 Porticos en sentido X del Bloque Estructural 67 
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tO. IS +• 11. IS* * 11.1S— * 13.10— ia.i i 

Portico W 

Figura 8 Porticos en sentido Y del Bloque Estructural 67 

Se encuentra la matriz de rigidez asociada todos los grados de libertad, 
indicados en la figura 9 y se condensa a las coordenadas laterales, con lo que se 
encuentra la matriz de rigidez lateral del portico. Aguiar (2014). 

Por un procedimiento similar se halla la matriz de rigidez lateral de cada 
uno de los porticos que conforman la estructura. 

Antes de terminar con la explication del calculo de la matriz de rigidez 
lateral de los porticos, conviene mirar los grados de libertad del portico V, en la 
figura 9 y darse cuenta que existe un grado de libertad horizontal a nivel de pisos, 
es decir se considera que las vigas son axialmente rigidas. El vector de colocation 
de los aisladores de los ejes 32, 33 y 36 tiene 6 cantidades diferentes de cero, en 
cambio de los aisladores de los ejes 39, 41 y 43 solo tiene 3 cantidades diferentes 
de cero ya que se encuentra empotrado en su base. 

El analisis sismico se realiza con inercias gruesas en los elementos, ya 
que van a trabajar en el rango elastico. Solo los aisladores sismicos incursionan 
en el rango no lineal por lo que se decidio trabajar con la rigidez efectiva del 
modelo de tres fases de Me Vitty y Constantinou (2015), esto con relation a la 
rigidez horizontal, la rigidez vertical se halla en funcion de la geometria del aislador 
interior. Constantinou et al. (2011). 



80 



Roberto Aguiar, Valeria Andrade, Byron Guaygua y Pablo Piedra 




10.15 v 11.15 >■ 11.15 13.10 * 13.11 

Figura 9 Grados de libertad de portico del Eje V del Bloque Estructural 67. 



7. ANALISIS SISMICO EN COORDENADAS DE PISO 

Para el analisis sfsmico se consideran tres grados de libertad por planta: 
dos componentes de desplazamiento horizontal y una rotacion de piso. Estos 
grados de libertad se ubican en el centra de masas de cada una de las plantas 
para que la matriz de masas sea diagonal. 

La matriz de rigidez en coordenadas de piso K E se encuentra con la 
siguiente ecuacion. 



K E =Y.A^ t K^A^ ( 2 ) 

Donde A (i) es la matriz de compatibilidad del portico i, que relaciona las 
coordenadas laterales de los porticos con las coordenadas de piso; es la 
matriz de rigidez lateral del portico i. 

En la figura 7 se aprecia que se tienen porticos de 3, 4, 5 y 6 pisos. Se 
cuenta como un piso al lugar donde se hallan los aisladores. Luego se van a tener 
matrices de rigidez lateral de 3 por 3; 4 por 4; 5 por 5, y 6 por 6. Para poder sumar 
matrices (ecuacion 2) todas tienen que ser del mismo orden, en el presente caso 
de 6 por 6. Por lo que las matrices de orden inferior deben ser completadas con 
ceros hasta tener matrices de 6 por 6. Pero los ceros no van en cualquier parte 
sino que van ubicados en la fila y columna donde no se tiene un determinado piso. 
Lo propio se hace con la matriz de compatibilidad de deformaciones A, para que 
tenga 6 filas y 18 columnas. 

En el libro Analisis Dinamico de Estructuras, con CEINCI-LAB (Aguiar 
2012) esta descrito el Metodo Modal Espectral. En efecto se presenta el 



Aplicaciones del Metodo Espectral 



81 



desacoplamiento del sistema de ecuaciones diferenciales y se deducen las 
ecuaciones con las que se halla los desplazamientos y fuerzas en cada modo de 
vibracion. 

Tambien se indica los criterios de combinacion modal, que mas se 
emplean. En este estudio se trabajo con el criterio CQC, Combinacion Cuadratica 
Completa, propuesto por Chopra (2001). 

Finalmente se destaca que los desplazamientos laterales se hallan con 
coordenadas absolutas, con respecto al movimiento del suelo. 

7.1 Propiedades dinamicas para sismo DBE y MCE 

En la tabla 11 se indica el periodo y el valor de los 6 primeros modos de 
vibracion, para el sismo DBE y en la tabla 12 para el sismo MCE. La estructura 
tiene 18 grados de libertad, en consecuencia 18 modos de vibracion e igual 
numero de perfodos de vibracion pero para no alargar la exposicion se presentan 
los valores de los 6 primeros modos. 



Tabla 11 Propiedades dinamicas para sismo DBE 





Modo 1 


Modo 2 


Modo 3 


Modo 4 


Modo 5 


Modo 6 


Periodo 


3.49 


3.46 


3.19 


0.75 


0.69 


0.57 


MODOS 


0.0005 


-0.0011 


0.0006 


-0.0001 


0.0019 


-0.0008 




0.0101 


-0.0331 


-0.0139 


0.0011 


0.0512 


0.0132 




0.0109 


-0.0355 


-0.0145 


0.0005 


0.0295 


0.0034 




0.0116 


-0.0375 


-0.0149 


-0.0004 


-0.0037 


-0.0091 




0.015 


-0.0414 


-0.0008 


-0.0013 


-0.0394 


-0.0107 




0.0179 


-0.0443 


0.0116 


0.0002 


-0.0663 


0.0077 




0.0011 


0.0004 


-0.0002 


-0.0019 


0 


0.0003 




0.033 


0.0129 


-0.0082 


-0.0525 


-0.0011 


0.0106 




0.0364 


0.0133 


-0.0049 


-0.0308 


0.0005 


0.0011 




0.0389 


0.0141 


-0.0051 


0.0033 


0.0002 


-0.0039 




0.0407 


0.0147 


-0.005 


0.0404 


0.0001 


-0.0057 




0.0419 


0.0148 


-0.0036 


0.0672 


0.0003 


-0.0025 




0 


0 


0 


0 


0 


-0.0001 




0.0004 


-0.0003 


0.002 


-0.0004 


0.0009 


-0.0025 




0.0004 


-0.0003 


0.0021 


-0.0003 


0.0007 


-0.0013 




0.0005 


-0.0004 


0.0022 


-0.0001 


0.0003 


0.0005 




0.0005 


-0.0004 


0.0023 


0.0003 


0.0001 


0.0029 




0.0005 


-0.0004 


0.0023 


0.0006 


0.0001 


0.0049 
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Tabla 12 Propiedades dinamicas para sismo MCE 





Modo 1 


Modo 2 


Modo 3 


Modo 4 


Modo 5 


Modo 6 


Perfodo 


3.57 


3.53 


3.25 


0.75 


0.69 


0.58 




0.0005 


-0.0011 


0.0005 


-0.0001 


0.0018 


-0.0008 




0.0102 


-0.0333 


-0.0137 


0.001 


0.0512 


0.0133 




0.0109 


-0.0356 


-0.0143 


0.0005 


0.0294 


0.0034 




0.0115 


-0.0376 


-0.0147 


-0.0004 


-0.0039 


-0.0091 




0.0149 


-0.0413 


-0.0005 


-0.0013 


-0.0395 


-0.0107 




0.0178 


-0.0441 


0.0118 


0.0003 


-0.0663 


0.0077 




0.0011 


0.0004 


-0.0002 


-0.0018 


0 


0.0003 




0.0332 


0.0128 


-0.0083 


-0.0525 


-0.0012 


0.0106 


MODOS DE 


0.0365 


0.0132 


-0.0049 


-0.0306 


0.0005 


0.0011 


VIBRACION 


0.0389 


0.014 


-0.0051 


0.0035 


0.0002 


-0.0039 




0.0406 


0.0145 


-0.005 


0.0405 


0.0001 


-0.0057 




0.0418 


0.0146 


-0.0036 


0.0672 


0.0004 


-0.0025 




0 


0 


0 


0 


0 


-0.0001 




0.0004 


-0.0003 


0.002 


-0.0004 


0.0009 


-0.0025 




0.0004 


-0.0003 


0.0021 


-0.0003 


0.0007 


-0.0013 




0.0004 


-0.0003 


0.0022 


-0.0001 


0.0003 


0.0005 




0.0005 


-0.0004 


0.0023 


0.0003 


0.0001 


0.0029 




0.0005 


-0.0004 


0.0023 


0.0006 


0.0001 


0.0049 



7.2 Resultados del analisis sismico en sentido X 

En la figura 7 se indico los porticos en sentido X, corresponden a los que 
tienen menos longitud y se hallan en la direction del voladizo. Por lo tanto, se trata 
del sentido mas critico. 

En la tabla 13 se presenta el factor de participacion modal r en los 6 
primeros modos de vibracion y las respectivas aceleraciones espectrales A d para 
el sismo DBE 

Ahora en la tabla 14 se indican los desplazamientos q que se hallan luego 
de aplicar el criterio de combination modal CQC y las aceleraciones de piso A v . 
Se destaca que los desplazamientos de la estructura se hallan multiplicando los 
valores que reporta el analisis espectral por Rcp p (p e -, para el sistema de aislacion 
son los que reporta el modal espectral. 

Los valores de la tabla 15 son similares a los de la tabla 13 pero para el 
sismo MCE y los resultados de la tabla 16 son similares a los de la tabla 14. 
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Tabla 13 Factores de participation modal y aceleraciones espectrales para sismo 



DBE en sentid 



oX 





Modo 1 


Modo 2 


Modo 3 


Modo 4 


Modo 5 


Modo 6 


r 


7.4929 


22.7786 


5.3027 


0.0087 


1 .4834 


0.1454 


Ad 

(m/s 2 ) 


1 .2397 


1.2514 


1 .3565 


9.7892 


10.5642 


12.7414 



Tabla 14 Desplazamiento y aceleraciones de piso para sismo DBE en sentido X 





Subsuelo 


Aislacion 


Piso 1 


Piso 2 


Piso 3 


Piso 4 


<7 

(cm) 


0.73 


28.09 


30.17 


31.88 


31.62 


31.07 


F 

(T.) 


1.18 


25.74 


74.53 


130.23 


137.54 


92.01 


A p 

(m/s 2 ) 


10.62 


1.26 


1.13 


1.10 


1.34 


1.66 



Tabla 15 Factores de participation modal y aceleraciones espectrales para sismo 

MCE en sentido X 





Modo 1 


Modo 2 


Modo 3 


Modo 4 


Modo 5 


Modo 6 


r 


7.5123 


22.7778 


5.2960 


0.0071 


1.4181 


0.1397 


Ad 

(m/s 2 ) 


1.7158 


1.7317 


1 .8786 


13.6744 


14.7591 


17.816 



Tabla 16 Desplazamiento y aceleraciones de piso para sismo MCE en sentido X. 





Subsuelo 


Aislacion 


Piso 1 


Piso 2 


Piso 3 


Piso 4 


q 

(cm) 


1.01 


40.88 


43.48 


45.93 


45.53 


44.71 


F 

(T.) 


1.55 


41.53 


107.33 


180.32 


186.84 


124.03 


A p 

(m/s 2 ) 


14.89 


1.72 


1.55 


1.52 


1.84 


2.26 



7.3 Resultados del analisis sismico en sentido Y 

En el sentido Y, la estructura tiene mejor comportamiento que el sentido X, 
debido a que tiene mas ejes resistentes y no tiene voladizo en ese sentido. Los 
resultados para el sismo DBE se indican en las tablas 17 y 18; para el sismo MCE 
en las tablas 19 y 20. 
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Tabla 17 Factores de participacion modal y aceleraciones espectrales para sismo 



DBE en sentid 



o Y 





Modo 1 


Modo 2 


Modo 3 


Modo 4 


Modo 5 


Modo 6 


r 


22.8262 


8.3773 


3.3475 


1.7139 


0.0028 


0.1501 


Ad 

(m/s 2 ) 


1 .2397 


1.2516 


1 .3584 


9.7892 


10.5642 


12.741 


Tabla 18 Desplazamiento y aceleraciones de piso para sismo DBE en sentido Y 




Subsuelo 


Aislacion 


Piso 1 


Piso 2 


Piso 3 


Piso 4 


q 

(cm) 


0.86 


23.79 


27.05 


28.94 


30.38 


31.66 


F 

(T.) 


0.38 


99.62 


126.69 


135.22 


111.16 


66.03 


A p 

(m/s 2 ) 


10.84 


1.31 


1.16 


1.13 


1.34 


1.67 


Tabla 19 Factores de participacion modal y ac 

MCE en sentid 


eleraciones espectrales para sismo 
lo Y 




Modo 1 


Modo 2 


Modo 3 


Modo 4 


Modo 5 


Modo 6 


r 


22.8527 


8.3169 


3.3547 


1.6391 


0.0041 


0.1431 


Ad 

(m/s 2 ) 


1.7158 


1.7317 


1 .8786 


13.6744 


14.7591 


17.816 


Tabla 20 


Desplazamiento y aceleraciones de | 


piso para sismo MCE en sentido Y. 




Subsuelo 


Aislacion 


Piso 1 


Piso 2 


Piso 3 


Piso 4 


q 

(cm) 


1.266 


37.40 


40.92 


43.60 


45.54 


46.88 


F 

(T.) 


0.884 


162.59 


185.65 


185.65 


142.409 


77.311 


A p 

(m/s 2 ) 


15.21 


1.78 


1.60 


1.57 


1.85 


2.27 



8. ANALISIS SISMICO PLANO 

Lo ideal es realizar un analisis sismico espacial, con tres grados de libertad 
por planta para incorporar el efecto de torsion pero es muy complicado presentar 
el procedimiento de calculo en forma detallada, por lo que en este apartado se 
presenta un programa general con el cual se realiza el analisis espectral piano, 
utilizando la librerfa de programas del sistema de computacion CEINCI-LAB su 
lectura sirve para entender el procedimiento de calculo. 

Se analiza un portico en sentido X (el 36) y un portico en sentido Y (el V), los 
resultados son parecidos a los que se hallan del analisis espacial. 
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8.1 Analisis de Portico en sentido X 

Se realiza el analisis sismico del Portico 36, cuya geometria y cargas 
consideradas se indica en la figura 10. Se destaca que las cargas verticales fueron 
obtenidas considerando la contribution de la carga de la losa y de las vigas 
perpendiculares a la direccion del portico. De tal manera que estan ligeramente 
mayoradas pero se justifica ya que se realiza un analisis sismico piano. 



7 448 T/m 




Figura 10 Portico 36 y cargas verticales consideradas. 

• Programa en CEINCI-LAB 

1 £ Programa para realizar el analisis sismico piano del portico 36 

2 % por el Metodo Modal Espectral 

3 £ Dr. Roberto Aguiar 

4 % 17 de diciembre de 2015 

5 % 

6 - load KLX36; % 5e carga la matriz de rigidez lateral del portico con ai3ladores 

7 - K=KLX36; 

8 - ml=0 .1*24.530; m2 =2 4.530; m3=2 7 . 705 ; m4=2 7.705; m5=2 7.705; m6=2 1.230; 

9 - M=[ml m2 m3 rn4 m5 m6 ] ; M^diag (M) 

10 - [T r f i r OM]=orden_eig (K r M) % periodos y modos de vibracion 

11 - zeda=0.2023; ^Factor de amortiguamiento de los aisladores, promedio 

12 - E= ( zeda/0 . 05 ) A Q . 3 ; % Para modificar el espectro en rango de aisladores 

13 - R=l; fip=l . 0; fie=0 . 9; IPara modificar espectro en rango estructura 

14 - to=2 . 0 ; ^Periodo a partir del cual se reduce espectro por B 

15 - [Ad] =espectro_necl5_mod_DE_l (R, f ip f f ie f B f to, T) ; ^Aceleraciones espectrale3 

16 - b=[l;l;l;l;l;l];% Factor de incidencia 

17 - Rfor i=l : 6 % Factores de participacion modal 

18 - gama { : , i> =abs ( ( ( f i ( : , i ) ) ■ *M*b > / { { f i ( : F i> ) ' *M* ( f i ( : r i } ) } ) ; 

19 - L end 

20 - na=6; % numero de grados de liber tad 

21 - [qt]=desplazamientos_modales_CQC [T r f i F Ad f gama,na,OM r zeda) %de3plazamientos 

22 - qine=R*f ip*'f ie*qt % Desplazamdentos Inelasticos de la estructura 

23 - [Yt, Ft ]=fuerza3_modale3_CQC (M F fi F Ad r gama,na,OM r zedaj % Fuerzas laterale3 

24 - [Fta r Ap]=fuerza3_modale3_aceleracion{M r f i r Ad r gama,na> % Aceleracione3 de piso 
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CEINCI-LAB es un sistema de computacion desarrollado por el primero de 
los autores de este articulo desde el 2008 y permite resolver las estructuras en 
forma amigable, para ello el usuario debe saber como se hace el calculo paso a 
paso; debe conocer la librerfa de programas que existen y saber los datos de 
entrada y los resultados que se esperan con cada programa. 

Pero no todo hace CEINCI-LAB, se deja a que el usuario programe algo 
del calculo, en el programa presentado se programo la determinacion de los 
factores de participation modal. 

En base a las cargas verticales de la figura 10, se obtuvo los pesos en 
cada piso y con esos datos se hallo la matriz de masas M, que aparece al principio 
del programa. La matriz de rigidez lateral se obtuvo con otro programa y aquf solo 
se carga el archivo que contiene ese matriz. 

Con CEINCI-LAB el usuario sabe a ciencia cierta que esta haciendo y 
tiene un mayor conocimiento del comportamiento de la estructura. 

El perfodo fundamental de vibracion que se encuentra es T = 3.57 s, para 
el sismo DBE, ligeramente mayor al que se hallo en el analisis espacial de T = 
3.49 s, esto se debe a que en el analisis piano se incremento la carga vertical, para 
ser un poco mas conservador el analisis. 

Los desplazamientos laterales, fuerzas laterales y aceleraciones de piso 
se presentan en la tabla 21 . 



Tabla 21 Desplazamientos, fuerzas laterales y aceleraciones de piso en 
Portico 36, ante sismo DBE 





Subsuelo 


Aislacion 


Piso 1 


Piso 2 


Piso 3 


Piso 4 


q 

(cm) 


1.09 


33.93 


33.31 


35.81 


37.56 


38.52 


F 

(T.) 


2.26 


25.57 


27.01 


29.36 


37.16 


36.50 


-4p 

(m/s 2 ) 


10.90 


1.51 


1.33 


1.26 


1.43 


1.71 



No se puede comparar directamente los valores de la tabla 21 con los 
valores de la tabla 14 debido a que en la tabla 21 se tienen desplazamientos en el 
portico p, mientras que en la tabla 14 se hallan desplazamientos en coordenadas 
de piso q. La relacion entre estos dos vectores viene dada por la matriz de 
compatibilidad de deformaciones A. 



P$ = A (n 4 (0 (3) 

Donde el fndice i hace relacion al modo de vibracion y el subindice j al 
numero de portico. De tal manera que para hallar los desplazamientos laterales en 
cada portico a partir de los desplazamientos y giros de coordenadas de piso, se 
debe encontrar la contribution en cada modo de vibracion y luego aplicar el criterio 
de combination CQC. 
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De todas maneras son comparables los valores de la tabla 21 con los de la 
tabla 14. Para encontrar las fuerzas laterales en cada portico, se hallo primero las 
fuerzas P® en cac | a modo con la siguiente ecuacion. 

P® = Ad ® r® M 0® 

Donde 4d® es la aceleracion espectral en el modo i; r® es el factor de 
participacion modal en el modo i; M es la matriz de masas; $® es el modo de 
vibracion i. 

Con las fuerzas laterales en cada modo de vibracion se encontro los 
cortantes de piso en cada modo V®, se aplico el criterio de combinacion CQC en 
los cortantes y a partir de estos cortantes se determino las fuerzas laterales por un 
procedimiento inverso a la forma de calculo de los cortantes. (Aguiar, 2012). Por 
esto el programa fuerzas_modales_CQC reporta los cortantes V t y las fuerzas F t . 

Las aceleraciones de piso, es un parametro fundamental para determinar 
las perdidas esperadas por cafda o volcamiento de los elementos que contiene la 
estructura para su determinacion en cada modo de vibracion se resuelve un 
sistema de ecuaciones lineales. Filiatrault (2015). 

M Af = P® (5) 

Donde Ap L) es la aceleracion espectral en el modo i. Luego se aplica un 
criterio de combinacion modal y se halla las aceleraciones de piso. 

8.2 Analisis de Portico en sentido Y 

En el sentido transversal se analizo el Bloque Estructural V, que se indica 
en la figura 1 1 . 
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El programa de calculo es muy similar al ya indicado, por lo que se omite 
su presentation. Los resultados obtenidos de desplazamientos horizontales, 
fuerzas laterales y aceleraciones de piso se indican en la tabla 22. 



Tabla 22 Desplazamientos, fuerzas laterales y aceleraciones de piso en 
Portico V, ante sismo DBE 





Subsuelo 


Aislacion 


Piso 1 


Piso 2 


Piso 3 


Piso 4 


<7 

(cm) 


1.12 


30.73 


30.26 


32.57 


34.20 


35.13 


F 

(T.) 


0.41 


38.88 


38.63 


45.35 


46.04 


31.83 


A p 

(m/s 2 ) 


10.06 


1.58 


1.42 


1.43 


1.68 


2.05 



Los resultados hallados son comparables con los indicados en la tabla 18. 
Destacando una vez mas que estos son a nivel de coordenadas de piso y los de la 
tabla 18 son a nivel del centra de masa. 

9. COMPARACION CON ESTRUCTURA SIN AISLACION 

Se realizo el analisis sismico espacial, en coordenadas de piso, del Bloque 
Estructural 67, suponiendo que no se va a colocar aisladores y los resultados a 
nivel de desplazamientos se compararon con los hallados en la estructura con 
aisladores sfsmicos, tanto para el sismo DBE como para el sismo MCE. 

En la figura 12 se aprecian los resultados para cuando el sismo actua en 
sentido X; la Ifnea con rojo representa el comportamiento de la estructura sin 
aisladores y la azul con aisladores. 

En la estructura con aisladores el mayor desplazamiento se tiene en los 
aisladores, como tal y la estructura en si no experimenta desplazamientos, con 
respecto al desplazamiento del aislador, menores a 5 cm., para el sismo DBE y 
menores a 8 cm, para el sismo MCE. 

En cambio en la estructura sin aisladores se tienen grandes deformaciones, 
alrededor de los 27 cm, para el sismo DBE y 38 cm, para el sismo MCE; esto 
implica dano en la estructura. 

En la figura 13, se compara el desempeno de la estructura en sentido Y, y los 
comentarios son similares al anterior. 

En la figuras 14 y 15 se presenta la deriva de piso, evaluada en el Centro de 
Masa, sin aisladores y con aisladores. En la curva correspondiente a la estructura 
con aisladores, existe una discontinuidad en el punto donde van los aisladores, ya 
que no tiene sentido dibujar la deriva de piso en los aisladores; en el subsuelo la 
deriva es muy baja y en los pisos superiores se incrementa en los pisos sobre el 
aislador pero en los ultimos pisos practicamente es nula. 
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Desplazamientos para sismo MCE en sentido X 




Figura 12 Comparacion de desplazamientos laterales en Centro de Masa, para sismo 
actuando en sentido X. Considerando estructura sin y con aisladores sfsmicos. 



Desplazamientos para sismo DBE en sentido Y 



30 




Desplazamiento (cm) 




Figura 13 Comparacion de desplazamientos laterales en Centro de Masa, para sismo 
actuando en sentido Y. Considerando estructura sin y con aisladores sfsmicos. 



Derlvas para sismo DBE en sentido X 




Oeriva de pi so en % 



Derivas para sismo MCE en sentido X 

29,64 
24,7 
19,76 

Can aislacion 

Sin aislacion 9 ' 33 

4,94 
O 




O 0,5 1 1,5 2 2,5 

Deriva de pi so en % 



Figura 14 Derivas de Piso en Centro de Masa, para sismo actuando en sentido X. 
Considerando estructura sin y con aisladores sfsmicos. 
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Derivas para sismo DBE en sentido Y Derivas P ara sism ° MCE en sentido Y 




Deriva de pi so en % Deriva depiso en« 



Figura 15 Derivas de Piso en Centro de Masa, para sismo actuando en sentido Y. 

Considerando estructura sin y con aisladores sismicos. 

En la estructura sin aisladores, para el sismo DBE la deriva de piso 
maxima es 1.6% y para el sismo MCE es 2.25%. Valores altos que implican dano 
en la estructura. El NEC 15 permite una deriva maxima del 2% para el sismo DBE, 
de tal manera que la estructura sin aisladores cumple con esta condition pero 
tener una deriva de piso de 1.6% implica dano en los elementos estructurales, no 
estructurales y sus componentes. 

En la estructura con aisladores sismicos, es 0.65% para el sismo DBE y 
0.71% para el sismo MCE, cantidades bajas. Por lo tanto, la estructura tendra muy 
buen comportamiento sismico. 

Retornando a los desplazamientos laterales, se debe manifestar que el 
desplazamiento lateral maxi mo encontrado para el sismo MCE es menor a los 48 
cm. En la construction se tiene un GAP de 60 cm. 



10.COMENTARIOS Y CONCLUSIONES 

Desde el punto de vista estructural y arquitectonico es muy buena option 
haber unido los bloques 6 y 7 de la UFA-ESPE. Con ello no se deja junta de 
construccion entre los dos bloques estructurales con lo que se aporta al proyecto 
arquitectonico y se compran 4 aisladores menos. 

Por otra parte, se realizo el analisis sismico en forma espacial y plana, 
empleando el Metodo Espectral. El comportamiento no lineal de los aisladores se 
obtuvo empleando el modelos de tres fases de Me Vitty y Constantinou (2015) que 
representa el estado del arte para los aisladores de triple pendulo de friccion. 

El analisis sismico piano, se lo realizo para poder mostrar el marco teorico 
y el uso del sistema de computacion CEINCI-LAB. 

Finalmente se comparo el comportamiento de la estructura con aisladores 
y de la estructura sin aisladores se observo que en la estructura con aisladores el 
desplazamiento lateral de los pisos es minimo, practicamente en la estructura 
aislada la superestructura se mueve como cuerpo rigido sobre los aisladores y 



Aplicaciones del Metodo Espectral 



91 



esto se refleja en las derivas de piso que son menores o igual al 0.65% para el 
sismo DBE y 0.71 % para el sismo MCE. 
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